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ABSTRACT 

We present new, high-to-intermediate spectral resolution stellar population models, 
based on four popular libraries of empirical stellar spectra, namely Pickles, ELODIE, 
STELIB and MILES. These new models are the same as our previous models, but 
with higher resolution and based on empirical stellar spectra, while keeping other in- 
gredients the same including the stellar energetics, the atmospheric parameters and 
the treatment of the Thermally-Pulsating Asymptotic Giant Branch and the Horizon- 
tal Branch morphology. We further compute very high resolution (R=20,000) models 
based on the theoretical stellar library MARCS which extends to the near-infrared. 
We therefore provide merged high resolution stellar population models, extending from 
~ 1000 A to 25,000 A, using our previously published high resolution theoretical 
models which extended to the ultraviolet. We compare how these libraries perform 
in stellar population models and highlight spectral regions where discrepancies are 
found. We confirm our previous findings that the flux around the V band is lower 
(in a normalised sense) in models based on empirical libraries than in those based on 
the BaSeL-Kurucz library, which results in a bluer B — V colour. Most noticeably the 
theoretical library MARCS gives results fully consistent with the empirical libraries. 
This same effect is also found in other models using MILES, namely Vazdekis et al. 
and Conroy & Gunn, even though the latter authors reach the opposite conclusion. 
The bluer predicted B — V colour (by 0.05 magnitudes in our models) is in better 
agreement with both the colours of Luminous Red Galaxies and globular cluster data. 
We test the models on their ability to reproduce, through full spectral fitting, the 
ages and metallicities of galactic globular clusters as derived from CMD fitting and 
find overall good agreement. We also discuss extensively the Lick indices calculated 
directly on the integrated MILES-based SEDs and compare them with element ratio 
sensitive index models. We find a good agreement between the two models, if the 
metallicity dependent chemical pattern of the Milky Way stars is taken into account 
in this comparison. As a consequence, the ages and metallicities of galactic globular 
clusters are not well reproduced when one uses straight the MILES-based indices, 
because subtle chemical effects on individual lines dominate the age derivation. The 
best agreement with the ages of the calibrating GCs is found with either element-ratio 
sensitive absorption-line models or with the full SED fitting, for which no particular 
weight is given to selected lines. 

Key words: galaxies: stellar content - stars: evolution - stars: fundamental param- 
eters. 



1 INTRODUCTION 

Galaxy evolution studies are reaching a high level of sophis- 
tication due to the very high quality of observational data 
permitted by modern technology, and the level of spectral 
details that such observations carry in (just to name a few 
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examples, the Sloan Digital Sky Survey - SDSS - Eisenstein 
al. 2011, Aihara et al. 2011; the GAlaxy Mass Assembly - 
GAMA - Driver et al. 2011; RESOLVE, Kannappan et al. 
2011; GMASS, Cimatti et al. 2008; the SINS survey, Forster- 
Schreiber et al. 2011). In this era of precision astrophysics, 
interpretative models, such as stellar population and galaxy 
models, need to keep pace with the fast observational devel- 
opment. 

Stellar population models allow the calculation of the 
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integrated spectral energy distribution (SED) of a synthetic 
stellar population - a star cluster or a galaxy - for arbi- 
trary choices of parameters such as the age, the Initial Mass 
Function (IMF), the metallicity, the star formation history, 
the Horizontal Branch morphology, etc. Models are based 
on three ingredients, namely the stellar energetics, the stel- 
lar spectra - onto which the energetics are distributed - and 
the computational method adopted for calculating the inte- 
grated model. At given energetics, a model SED of higher 
spectral resolution is obtained by employing higher resolu- 
tion stellar spectra. In this paper we present the higher- 
resolution version of the Maraston (2005) stellar population 
models. 

In the Maraston (2005) models, the library of 
low-resolution , theo r etical stellar spectra compiled by 
iLeieune etaU l| 19971 . | 1998T ). was adopted. This is l argely 
based on the classical model atmospheres of iKuruczl (|l979T ) 
and a patchwork of different calculations for cool, M-type gi- 
ants and dwarfs. The final assembly was further corrected in 
order to match empirical relations between stellar temper- 
ature and colours. The library does not contain Thermally- 
Pulsating Asymptotic Giant Branch (TP-AGB) carbon-rich 
and oxygen-rich stars, for which in Maraston (2005) the em- 
pirical, time-average d spec tra of C- and O-type stars from 
lLancon fc Mouhcind <|2002T l were utilised. 

A higher spectral resolution is required for a detailed 
modelling of emission and absorption lines. This is not the 
only necessary feature because element abundance ratios 
play an equally important role, but it is the first obvious step 
to resolve the lines properly. In turn, a high spectral resolu- 
tion SED can be obtained either with theoretical or empiri- 
cal stellar spectra. The two approaches have both pros and 
cons. While theoretical spectra give in principle high flexibil- 
ity on the stellar parameter coverage (temperature, gravity, 
metallicity) and in the wavelength extension, known prob- 
lems exist in th e modelling of spec if ic lines (e.g. the Hg line 
at ~ 4861 A, ICoelho et all 120071 ; IWalcher et ' all l2009h as 
well as con tinuum regions due to line blanketing (around 
the V band.lWorthevlll994l ) and around the g, r SDSS bands 
|Maraston et al.ll2009al . see below) . Empirical stellar spectra 
have all lines that are observed in stars, but their intensity 
and flux ratios depend on the chemical enrichment history of 
the site where the star is found. Moreover, empirical libraries 
usually cover a limited region of stellar parameter space and 
a limited A range. Hence, the approach of combining empir- 
ical and theoretical spectra is the most convenient one. 

The first version of high-resolution Maraston-type mod- 
els was presented in iMaraston et al.l l|2009bl ), where fully 
theoretical, very high-resolution SEDs in the mid- and far- 
UV (A < 4350) were published. At optical wavelengths, sev- 
eral empirical stellar libraries have been made available in 
the recent years, containing enough stars for nearly-complete 
stellar population models in a relatively broad wavelength 
range to be calculated (the Pickles library, Pickles 1998; the 
Jones's library, Jones 1999; the STELIB library, Le Borgne 
et al. 2003; the ELODIE library Prugniel & Soubiran 2001, 
and the MILES library, Sanchez-Blazques et al. 2006) . Stel- 
lar population models based on these libraries exist in the 
literature (e.g. Bruzual & Chariot 2003 for the STELIB li- 
brary; Vazdekis et al. 2010 for the MILES library). 

Our preliminary attempt with these libraries 
l|Maraston et al.l l2009al ) allowed us to improve upon a 



long-standing mismatch between the colours of luminous 
red galaxies in SDSS and stellar p opulation models based 
on theoretical stellar spectra l|Eisenstein et all 1200 ll ; 
I Wake et alj|2006j ). Also, a few remaining problems with the 
broadband colours of stellar population models, such as 
the (B-V) colo urs of Milky Way GCs (e.g. IWorthevi Il99l : 
IMarastonl 2005) could be connected to the same origin, and 
we shall show in this article that this was indeed the case. 

In this paper we present the full grid of M05-type mod- 
els for a number of these different optical empirical libraries. 
In addition, we also employ the theoretical MARCS library 
for modelling the spectra of old populations with a high res- 
olution up to the near-IR. We shall test the models with GCs 
data, and compare them with results from other authors. 

The paper is organised as follows. In Section [2] we re- 
call the main features of the evolutionary population syn- 
thesis code and in Section [3] those of the empirical spectral 
libraries. In Section|4]we describe the implementation of the 
empirical stellar spectra in the synthesis code, and also the 
calculation of models based on high-resolution theoretical 
spectra. Stellar population models are presented and amply 
discussed in Section [5] including a comparison with the lit- 
erature, while in Section [6] the model fit to globular cluster 
data is treated. A summary is given in Section [7] 



2 SUMMARY OF THE POPULATION 
SYNTHESIS CODE 

The foundation upon which these mod e ls have been built 
was thoroughly described in IMarastonl (|l998l . 120051 , here- 
after M98 and M05, respectively). Here, we present a brief 
recapitulation of those models. 

Evolutionary population synthesis models assume a 
stellar evolution prescription, which i n this case consisted 
of the isochrones and stellar tracks bv lCassisi et al.l (Il997l) 
for ag es larger than ~ 30 Myr and by ISchaller et al.l 
(1992) for younger populations. Sets of models have been 



also computed with the Padova stellar evolutionary models 



dGirardi et al. |2000D. The th eoretical fuel for the TP-AGB 
phase - from Renzinil |l992) - was calibrated with obser- 
vational data. Reimers-type mass-loss was applied to the 
Red Giant Branch track in order to generate blue and inter- 
mediate Horizontal Branch (HB) morphologies as observed 
in Milky Way globular clusters, and also blue HBs at high 
metallicity. 

The s ynthesis techn ique exploits the fuel consumption 
approach (jRenzinil [1981) for calculating the energetics of 
post Main Sequence phases, i.e. the amount of fuel avail- 
able for nuclear burning using the evolutionary track of the 
turnoff mass. In the M05 models, the fuel is used as integra- 
tion variable in post Main Sequence. For the Main Sequence, 
the standard isochrone synthesis technique is applied, where 
the integration variable is the mass. 

As for the st ellar s pectra, the BaSeL assembly 
l|Leieune et al.lll997h of the iKurucj (|l979l ) library was em- 
ployed for all evolutionary phases except the TP-AGB, for 
which empiric al carbon-rich and oxygen -rich stellar spec- 
tra were used (lLancon fc Mouhcind l2002h . Finally, the code 
structure is modular (see Maraston 1998), i.e. the three in- 
gredients - the energetics, the atmospheric parameters tem- 
perature and gravity, and the individual stellar spectral - are 
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stored in independent matrices, which allows one to easily 
experiment on their differential impact on the final mod- 
els. This structure was particularly useful to compute the 
models of this paper for which we only vary the spectral 
transformation matrix. 



3 THE LIBRARIES OF EMPIRICAL STELLAR 
SPECTRA 

Four different libraries of flux-calibrated empirical ste l- 
lar spectra have been considered, namely iPicklesl |l998), 
ELODIE, STELIB, and MILES. A brief summary of the 
properties of these libraries can be found in Table Q] Note 
that the values of T e //, log(g), and [Fe/H] reported there are 
minimum and maximum values only, and therefore say little 
about the actual sampling of the parameter space. For exam- 
ple, in the MILES library only three stars have a metallicity 
higher than [Fe/H]= +1.0. The actual range interesting for 
population synthesis modelling is illustrated in Figures [6][9] 
In Table [T] we provide two values for the spectral resolu- 
tion, the nominal one cited in the reference papers and the 
actual one as was determined by later works. For MILES wc 
found that the actual resolution of the library is somewhat 
coarse r tha n the 2.3 A st a ted in ISanchez-Blazauez et al.l 
20061 and IVazdekis et all l|20ld ) and is instead 2.54 A 
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_ 201 lh . This res ult has been independentl y 
confirmed by the MILES group ( Falcon-Barroso et al.ll201l[ ) . 
For STEL IB. iMacArthur et al.l ll2009h found that the resolu- 
tion of the lBruzual fc Chariot! \200$ ) models is coarser than 
the resolution of individual STELIB spectra, which is the li- 
brary adopted in those models (see Section 3.2). The actual 
resolutions should be used. 

Below follows a review of important features of each 
library, together with a description of any alterations we 
may have imposed upon the spectra. 

Note that stars from the galactic Bulge are not included 
in these libraries, and also not spectra of carbon- and/or 
oxygen-rich TP-AGB stars. For the latter, we use the same 
empirical spectra as in the M05 models (see Section 3.6). 

3.1 The Pickles library 

The Pickles library (|Picklesl Il998h differs from the other 
three in that it does not contain spectra of individual stars. 
Rather, the author has vacuumed the literature for a versa- 
tile collection of stellar libraries and combined the various 
observations, in order to provide standard spectra for all nor- 
mal spectral types and luminosity classes. Because they are 
(for the most part) averages of spectra from several sources, 
the Pickles spectra should not be as affected by systematic 
and random observational errors, as a library containing sin- 
gle observations of each star using the same setup. 

The library has exemplary coverage of evolutionary 
phases, see Fig[l] from the coolest M dwarfs and giants 
(T e ff ~ 2500 K) to hot O-type stars and supergiants. Only 
solar metallicity models can be computed from this library, 
however, since there are not enough spectra to cover all evo- 
lutionary phases to a satisfactory degree for non-solar metal- 
licities. The spectral resolution is R= 500 and the wave- 
length range is very wide, ranging from the far-UV (1150 A) 
to the K-band (2.5 urn). The extension into the UV in the 
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Figure 1. Evolutionary coverage at Zq for the Pickles library. 
Superimposed are the 1 M yr, 1 Myr and 12 Gyr is ochrones from 
iLeieune fc Schaerer (2001) and Girardi et al. (20q3), respectively. 



Pickles standard spectra was accomplished using the low- 
dispersion Intern ational Ultraviolet Explorer (IUE) spectra 
jHeck et alJll984h . Above ~1 urn, about half of the spectra 
lack spectroscopic observations, and the author constructed 
a smooth energy distribution from broadband photometry 
here. This may imply that near-IR absorption features may 
not be well resolved. Also note that, at these wavelengths, 
the M giant spectra are synthetic. 

The effective temperatures of the Pickles spectra were 
assigned by means of ( V-Ic) and ( V -K) colour- temperature 
relatio ns mainly from iBesselll (|l979h , and from Fluks et al.l 
jl994 ) for the M giants. No accurate metallicity calibration 
has been performed, other than sorting the spectra into so- 
lar or non-solar metallic i ty ba sed on a recipe described in 
iPickles fc van der Kruitl (|l990l 'l. 



3.2 STELIB 



iLe Borgne etTafl |2003) obtained roughly 250 stellar spec- 
tra of fairly high resolution (~ 3 A FWHM) and with good 
wavelength coverage. Most spectral types and luminosity 
classes are included, with a fair coverage in metallicity. There 
is a deficit of stars with [Fe/H] < —1.0. Several spectra 
lack atmospheric parameter determinations, and many do 
not have complete spectroscopic observations over the en- 
tire quoted wavelength range. In addition to this, all spec- 
tra available to the public suffer from telluric contamination. 
We have circumvented this problem by simply interpolating 
over the most affected wavelength regions, more precisely 
6850-6950 A, 7580-7700 A, and 8850-9050 A. Granted, this 
might not be the most aesthetically pleasing solution, and 
it is by no means a complete eradication of all telluric ab- 
sorption features, but it acts as a warning beacon to users 
and ensures that no science can be done within said regions. 
The STELIB library has previously been adopted for stellar 
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Table 1. Properties of the empirical libraries 



Reference 


No. of 


Wavelength 


nom. FWHM 


actual FWHM 


T e// (K) 


log(g) 


[Fe/H] 




stars 


range (A) 


(at 5500 A) 


(at 5500 A) 








Pickles (1998) 


131 


1150-25000 


11 A 


11 A 


2500-39800 


[-0.5, +5.5^ 


[-0.8, +0.6 ] 


STELIB 


249 


3200-9300 


3.0 A 


3.4 A 


3250-51000 


[0.0, +5.0] 


[-2.75, +0.6] 


MILES 


985 


3500-7430 


2.3 A 


2.54 A 


2750-36000 


[0.0, +5.5] 


[-2.9, +1.6] 


ELODIE v3.1 


1388 


3900-6800 


0.55 A 


0.55 A 


3100-51000 


[-0.25, +4.9] 


[-3.0, +1.0] 



a Inferred from matching to the isochrones of lGirardi et al ] [|2000l l in a My-(B-V) CMD. 



population synthesis bv lBruzual fc Charlotj (|200l hereafter 
BC03), to which we refer for more information. These au- 
thors replaced the contaminated segments with theoretical 
Kurucz spectra for stars cooler than 7000 K, and with Pick- 
les spectra for hotter stars. 

Upon inspec tion of the indi v idual spectra of the 
STELIB library, iMacArthur et all (|2009h found an rms 
spread of 1.09 A in starting wavelength and 0.00033 A in 
wavelength dispersion. If not corrected for, this will cause 
the coaddition of any spectra into the SED of a stellar pop- 
ulation to have a degraded resolution compared to the reso- 
lution of the individual stellar spectra, acting effectively as 
a velocity broadening. In addition to this, they also found 
an rms error of 11.4 km s - 1 for the radial velocities adopted 
bv lLe Borgne et al l (j2003T l. and small systematic deviations 
in the wavelength scale, of no more than ±0.5 A. Both of 
these attributes will further add to the downgrading of the 
resolution of the stellar population model. 

In order to rectify these problems with the STELIB 
wavelengths, thus avoiding loosing valuable resolution, we 
have applied small, non-linear shifts to the wavelength scale, 
derived for each individual STELIB spectrum, by fixing the 
positions of eight strong absorption features, and use inter- 
polation to obtain the wavelength shifts for all flux points 
in between. The spectra were then rebinned to our adopted 
scale of 3201-9296.5 A in steps of 0.5 A. The reason for 
straying from the 1 A step-size is that we found that the 
rebinning was accompanied by a small, but non-negligible, 
smoothing of some absorption features, which could be made 
insignificant with a finer sampling. 

The majority of atmospheric parameter s in STELIB 
were taken from the ICavrel de Strobel et all (|l997l , 120011 1 
catalogues of [Fe/H] determinations, which are compiled 
from a vast number of different sources in the litera- 
ture. These values have be en complemented with additiona l 
T e ff determination s fro m iBT ackwcll & Lvnas-Gravl (|l998h . 
Idi Benedettol (|l998T l and lAlonso et alj I l|l996l . ll999h . Colour- 
temperature relations from the latter have also been used 
where available. Several stars, that are a lso part of the 
ELODIE library (jPrugniel fc Soubiranl200lT) . were given the 
atmospheric parameters listed therein, if no other values 
could be found. If more than one determination of a parame- 
ter for the same star was available, these were averaged, with 
a larger weight given to the most recent determinations. 

3.3 MILES 



The MILES library (|Sanchez-Blazquez et"afll2006to provides 
spectra for almost 1000 stars, covers most evolutionary 



stages - based upon what can be expected to exist in our 
Galaxy - at all metallicities down to [Fe/H] ~ —2.0. Only 
stars at the extreme ends of the temperature scale are lack- 
ing, such as giants cooler than around 3400 K, and young 
stars hotter than ~26000 K. The resolution is high (2.54 A 
FWHM), and the wavelength range is primarily confined to 
the optical (3500-7400 A). Values for all atmospheric param- 
eters are taken from a large number of sources in the liter- 
ature, but a careful attempt has been made to homogenise 
the different scales. The reader is referred to lCenarro et al] 
(2007) for a detailed description of this procedure. When 
it comes to the effective temperatures, however, it should 
be noted that the reference system adopted for this cali- 
br ation comprises only values within the range 4000-6300 
K l|Cenarro et al.ll2007l L This may imply that the temper- 
atures of stars outside this approximate interval may not 
be as homogeneous, even if an effort is made at tying every 
temperature to the system defined in the above interval. 

Linking to t he discussion in the Introduction, 
iMilone et al.l (|2010t ) published [Mg/Fe] abundance ratios for 
roughly 75 % of the stars in the MILES library. The char- 
acteristic 'knee' (see their Figure 1), which signals the onset 
of enrichment from SN la, occurs at [Fe/H]~ —0.9, below 
which the [Mg/Fe]-enhancement lies fairly constant around 
+0.4, albeit with a large scatter. Stellar population models 
incorporating MILES spectra will therefore follow a similar 
abundance pattern. The same will happen with the other 
libraries in the metal-poor regime. 

3.4 ELODIE 

Co mpared to the original libra ry of ELODIE echelle spec- 
tra dPrugniel fc Soubiranl 2001), the latest version, ELODIE 
v3.1 (|Prugniel et alj|2007i ). contains roughly twice as many 
stars and a slightly increased wavelength range. The latter 
is still rather narrow, spanning only around 3000 A in the 
optical, which is not surprising considering the very high 
resolution (R~10,000). One may argue about the relevance 
of applying stellar population models of such high resolu- 
tion to velocity-broadened galaxy spectra, and for this pur- 
pose alternative versions of the library with lower resolutions 
may be computed. There are, however, stellar populations in 
which velocity broadening is less prominent, such as dwarf 
galaxies and star clusters, that might benefit from such a 
high-resolution model. 

Most stars in the ELODIE library have solar or near- 
solar metallicity. Only 68 stars out of 1388 are listed as 
having [Fe/H]< —1.0, which limits our ability to pro- 
duce ELODIE-based models at low metallicity. In the origi- 
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nal paper l|Prugniel fc SoubiranllioOll) at mosphe ric parame- 



1998). or taken from the literature (ICarnev et al.l 


1994; 


Alonso et al.l 1996L 199S 


; Cavrel de Strobel et alj 


19971, 


200 ll; Thevenin 1 19981; 


Blackwell & Lvnas-Gravl 1 


1998TI 



1 19991 ) were also used. Final parameter values were taken 
as averages of multiple determination s, giving lower weigh t 
to old references. In the latest version i|Prugniel et alj|2007l ). 
these parameters have been updated using more recent pub- 
lications. 

The stellar parameter determinations and flux calibra- 
tion of the stars in the ELODIE library have also been as- 
signed quality flags by the authors. We take advantage of 
this to create as reliable SSP models as possible, by allow- 
ing only stars with the highest quality of flux calibration to 
contribute, and in cases where two or more stars occupy the 
same location in stellar parameter space, choosing the one 
with the most reliable parameter estimations. Some stars 
have also been observed multiple times. Again, in such cases 
we chose the spectrum with the highest values of the quality 
flags. If there are two or more spectra of the same star and 
with the same quality flags we compute an average spectrum 
for this star, taking also the average of the stellar parame- 
ters. 

3.5 The flux calibration 

The STELIB, MILES, and ELODIE libraries have a few 
stars in common, which gives us the possibility to obtain 
some indication of the quality of flux calibration. In FigfJ] 
we compare three stars of intermediate temperatures, chosen 
to represent a large range in metallicity and gravity. Spectra 
have been smoothed to a common, low resolution of R=500, 
and normalized to unity at 5050 A, in order to make any 
differences in the continuum slope more clearly distinguish- 
able. None of the stars are known to be a variable. As can 
be seen in the figure, for at least two of the stars there are 
surprisingly large discrepancies, hinting at the uncertainties 
involved in accurately flux calibrating a stellar spectrum. 
If the differences in flux calibration are due to systematics, 
these will also be transferred over to the SEDs of the stellar 
population models. Clearly visible are the telluric absorption 
features in the STELIB spectra around 6900 A. 

A careful description of the complex flux cali- 
bration process for the ELODIE spectra is given by 
IPrugniel fc Soubiranl (120011 ). Since this library is a compi- 
lation of archive spectra rather than being part of a spe- 
cific observation programme, no standard stars have been 
observed. The instrumental response is therefore evaluated 
by comparisons to several other low- and high-resolution li- 
braries, and the spectra are corrected thereafter. At 3900 A 
- the blue limit of the ELODIE spectrograph - the global 
efficiency of the instrument drops by a factor of 50, which 
may compromise an accurate flux calibration in the blue end 
of the spectra. 

The flux cal i bratio n of the STELIB stars was tested in 
iLe Borgne et all (|2003l ) by comparing synthetic photome- 
try against published p hotometry in the Lausanne database 
i|Mermilliod et al.lll997h . The (B-V) colour is by far the one 
with the smallest rms dispersion, 0.083 mag, with no dis- 
cernable offset. The (U-B) and (R-I) have the largest dis- 



persions, which is due to inhomogeneities in the photometric 
systems, and the fact that the STELIB stars do not have the 
required wavelength coverage for the U- and /-band filters. 
There appears to be a slight trend in the (V-R) colour, such 
that the reddest (coolest) stars appear somewhat bluer in 
the STELIB library. It is argued by the authors, that this 
is most likely due to differences in the photometric systems. 
It should also be noted that these tests were performed on 
stars not corrected for interstellar reddening. All this taken 
together makes it difficult to assess the accuracy of the flux 
calibration in bands other than B and V. 

The MILES stars have go ne through a rigorous 
spect rophotometric quality control (|Sanchez-Blazquez et al.l 
2006). An internal consistency check, comparing several ob- 
servations of the same stars, revealed that the random rms 
error affecting the flux calibration is around 0.013 mag 
for the (B-V) colour. An extern al comparison w i th th e 
Lausanne photometric database (Mcrmilli od et all Il997l ). 
showed a mean offset of 0.015 mag for (B-V), in the sense 
that the MILES spectra are bluer. The rms dispersion of 
0.025 mag is con s idera bly smaller than the one obtained in 
ILe Borgne et all l|2003l ) for the STELIB stars. For stars in 
common with the STELIB library, the MILES (B- V) colours 
are, on average, 0.01 mag bluer, but with a trend that es- 
sentially all stars with T e // <5600 K are bluer, whereas al- 
most all stars hotter than ~7200 K appear redder (see their 
Fig. 11). The authors also defined seven box filters with effec- 
tive wavelengths between 4000 and 6200 A, in order to com- 
pare synthetic colours for stars in common with the STELIB 
and ELODIE libraries. For all possible combinations of fil- 
ters allowed by the wavelength range in common between the 
libraries, any statistically significant offsets show the MILES 
stars to exhibit bluer colours. For colours involving the fil- 
ters with the largest central wavelengths, there once again 
appears to be a trend of the coolest stars being significantly 
bluer in the MILES library (their Fig. 9), whereas the sit- 
uation is the opposite for the hottest stars. It is difficult 
to assess the flux calibration red- wards of around 6500 A, 
since this portion of the spectra is more or less left out of 
the quality check. 

The trends summarised above will be mirrored in the 
SSP models (cf. Section 5). 



3.6 On the temperature scales of empirical 
libraries 

In this work we opt to retain the stellar parameters as 
they have been published by the authors of the various li- 
braries, because our aim is to assess the joint effect of a 
library /temperature scale/stellar parameters on stellar pop- 
ulation models, with the ultimate goal at understanding 
different results that are obtained on galaxies using differ- 
ent models. As discussed, the temperature scales are some- 
what inhomogeneous, which may cause differences in the 
integrated stellar population spectra (e.g. Percival & Salaris 
2009) that can sometimes be hard to explain. For example, 
we could attribute to the adopted temperature scale for red 
giants part of the reason why the near-IR flux of MILES 
stars drop longward 6000 A (see Appendix A). On the other 
hand, a dedicated publication is required in order to homo- 
geneise the different temperature scales. Star numbers vary 
a lot between different libraries, hence a common scale built 
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Figure 2. Comparison of the STELIB, MILES, and ELODIE libraries. The three panels depict the spectra of three 
stars present in all libraries, smoothed to a common, low resolution of R=500, where a) is a metal-poor sub-dwarf; 
b) a super-solar metallicity MS star; and c) a metal-poor giant 



avoiding extrapolation could be worked out only for a sub- 
sample of stars, with the effect of reducing the stars available 
for population synthesis. 

In order to get a sense for the effect of different temper- 
ature scales, we used MILES, that contains more stars than 
the other libraries, and placed the STELIB library on the 
MILES temperature scale. This was achieved by measuring 
the B — V colours of the MILES stars at a given metallic- 
ity and fit a relation between B — V and temperature. A 
MILES-based temperature scale for STELIB was then ob- 
tained by deriving the temperature of STELIB stars using 
their B — V . The resulting relation is shown in Figure [3] 
The effect on the stellar population model will be shown in 
Section [5^31 



a coarser resolution than most of the adopted libraries 
(R ~ 1100). For keeping the TP-AGB phase in the new 
models, we have re-binned the TP-AGB spectra to the 
sampling proper of the various empirical libraries. Except 
for the Pickles-based models, this implies that at the ages 
relevant to the TP-AGB phase in the M05 models (0.2 
to 2 Gyr) and at A > 6000 the actual resolution of the 
models will be slightly coarser. Note that - because we 
have not downgraded the TP-AGB spectra - there is an 
overall improvement in the definition of the carbon and 
oxygen-rich star spectral breaks in the integrated model 
compared to M05. 



3.7 Keeping the TP-AGB empirical spectra from 
M05 

The TP-AGB empir ical stellar spectra by 
lLancon fc Mouhcind |2002r i used in the M05 code have 



4 IMPLEMENTATION OF THE EMPIRICAL 
SPECTRA 

One of the main advantages of using theoretical spectra is 
that they can be computed to sample the stellar parameter 
space arbitrarily well. Inserting them into a stellar popu- 
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Figure 3. Fit to the B — V colour of MILES stars vs effective 
temperature. 

lation synthesis code is therefore relatively straightforward. 
When dealing with the more sparsely sampled grid of stel- 
lar spectra that are found in the empirical libraries, however, 
one has to use a slightly different approach. 

4.1 Assigning empirical stellar spectra to the 
stellar track 

In the M05 code, a stellar spectra is assigned to each [T e //, 
log(g)]-bin by means of a quadratic interpolation in T e // 
and log(g) at given metallicity. This approach is generally 
not feasible when empirical libraries are used, due to their 
coarser sampling of the parameter space. 

Instead, for each empirical library we created differ- 
ent tables of stellar spectra, corresponding to the differ- 
ent evolutionary phases, i.e. MS, SGB, RGB, HB, E-AGB, 
and supergiants. The categorisation was performed based 
mainly on the surface gravities and luminosity classes of the 
spectra, but other discriminants, such as visual inspection 
of the spectra, was used in uncertain cases. More specifi- 
cally, the cuts in surface gravity that were applied for each 
respective evolutionary phase are: MS (log(g)>4.0), SGB 
(3.5<log(g)<4.0), RGB, HB and E-AGB (1.0<log(g)<3.5), 
and supergiants (log(g)<1.0). These cuts are, of course, only 
approximate and depend somewhat on temperature as well; 
it is not unusual for hot MS stars (spectral class A or B) to 
have surface gravities below 4.0 dex. Spectra on the border- 
line between two evolutionary phases were therefore permit- 
ted to be present also in both corresponding tables. 

Within one evolutionary phase, the actual representa- 
tive spectrum for each mass bin with its given temperature 
is then calculated by means of linear interpolation of the log- 
arithmic fluxes in log(T e //). The number of stars assigned 
to each phase is never smaller than 10. We have checked the 
stability of the procedure by using fewer or more stars and 



fixed the procedure when the results were found to be stable. 
Note that the whole procedure is also checked at the bot- 
tom of the calculations by comparing the integrated models 
with other calculations based on theoretical stellar spectra 
or with literature calculations (see Section [5j. 

4.2 Scaling the empirical stellar spectra 

Another difference between theoretical and empirical stellar 
spectra is the flux units. Whereas theoretical spectra prefer- 
entially are given in absolute units, such as luminosity (ergs 
s _1 A -1 ) - or any other unit that can easily be converted 
into luminosity; the main point being it makes the popula- 
tion synthesis very straightforward - things are more com- 
plicated with empirical spectra. Either fluxes are given as 
they were measured (ergs s _1 cm" 1 A -1 ) and thus depen- 
dent on the distance to the stars, as is the case with STELIB 
and ELODIE spectra, or they have been normalised to unity 
at some wavelength, a method adopted in the Pickles and 
MILES libraries. In both cases, the true relative energy scale 
between the spectra is washed out. 

In order to be able to use a consistent approach for 
all the empirical libraries, the following prescription was 
adopted. 

i) We (re-) normalize all stellar spectra to the averag e 
flux in a 100 A passband around 5550 A, as in iPicklesI (|l998h . 
This choice is not crucial for the different libraries at differ- 
ent resolution because we use the average flux within the 
band. 

ii) The spectra are then scaled with the luminosities 
(again using the average value in a 100 A passband centered 
at 5550 A) of theoretical spectra with the same stellar at- 
mosph eric parameters (inte rpolated as necessary) from the 
grid of iLeieune et al.l (|l997T ) - which were the input of the 
M05 models. In short, the scaling of the empirical spectra is 
based on theoretical spectra, but our procedure allows us to 
keep, as much as possible, the continuum shape of the empir- 
ical spectra intact. For the MILES case we have also verified 
that the choice of width and wavelength is not crucial, by 
recalculating the scaling using a much larger A range, as 
adopted by other authors (see Section 5.3). 

We note already now that the adopted procedure de- 
termines important features of the integrated stellar popu- 
lation model, a point on which we return when comparing 
our models with others from the literature (Section 5.3). 

4.3 Topping up missing cool dwarfs 

Empirical libraries consistently lack cool dwarfs (precise 
temperatures naturally vary with metallicity, but rough 
guidelines are < 5000 K for STELIB, ~ 4000 - 4500 K for 
ELODIE, and ~ 4000 K for MILES at non-solar metallici- 
ties), hence a supplement of theoretical spectra is necessary 
for completing the whole isochrone. For our purposes, the 
latest version of the MARCS library l|Gustafsson et alfeOOSt ) 
was deemed the most suitable, comprising a comprehen- 
sive grid of cool stars (2500-8000 K) at very high resolution 
(R=20,000) over a large wavelength range (1300-200000 A). 
First, the standard chemical composition of these spectra 
follow the general trend of Milky Way stars as the empirical 
libraries, i.e. display a moderate enhancement of a-elements 
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([<*/Fe]= 0.1-0.3) at [Fe/H] between -0.25 and -1.0, with 
more metal-poor objects having [a/Fe]= +0.4.0 Quite rel- 
evant, we also verified that the optical colours of MARCS 
stars and the corresponding integrated SSP models based 
solely on MARCS spectra are in good agreement with those 
based on empirical stars only (Section 0. This implies that 
our addition of these cool dwarfs does not alter the empirical 
spectral shape (around the U-band). We have also verified 
that the strength of integrated Lick indices is not affected by 
the use of cool MARCS stars, in spite some lines do not agree 
with the correspondent features as obtained from the extrap- 
olations of MILES-based fitting functions (Johansson et al. 
2010) into the very low-temperature regime. This happens 
as very cool dwarfs and giants {T c g < 3500 °K) have little 
effect on integrated optical line-strengths, unless the IMF is 
bottom-heavy (Maraston et al. 2003). It is well known that 
spectral lines in empirical and theoretical libraries do not 
always agree (e.g. Tripicco & Bell 1995; Korn, Maraston & 
Thomas 2005), but a proper comparison at very low tem- 
peratures has not yet been done. Since - as we mentioned 
- most empirical libraries lack very cool stars, one cannot 
fully trust that the extrapolation of the lines observed in 
warmer stars is correct. It goes beyond the scope of this pa- 
per to examine each spectral feature of the MARCS spectra 
in detail, and we defer this to a future publication. 

Before entering the synthesis, the MARCS spectra were 
smoothed to the lower resolution of the empirical libraries. 

We have also calculated SSP models fully based on the 
MARCS library (Sec. US). 



4.4 Removing broken spectra 

During the calculation of the MILES-based models a com- 
parison to the value o f Hp from the fitting functions of 
Ijohansson et al.l (|201dh was performed, as we noted a dis- 
crepancy in the corresponding models. We found that in- 
deed two stars, that otherwise do not display any pecu- 
liarity regarding the shape of their SEDs, had conspicuous 
Hp strengths, one even had some Balmer emission. These 
were removed from the synthesis and replaced with stars 
of similar SED shape and normal Balmer strength values. 
The good agreement we find between solar-scaled models 
of Lick indices base d on index-by-index fitting functions 
|Thomas et al.ll201 if) and those calculated on the full SED 
presented in this paper comes from removing these anoma- 
lous stars (see Section 15.5 [I , This example supports our view- 
point that - when dealing with absorption lines using obser- 
vations - it appears to be safer to take a look at the in- 
dividual lines for spotting anomalous outliers, and the FF 
method is an effective way of doing so. 



4.5 £/V-extension of the models 

For galaxy formation and evolution studies, an exten- 
sion into the ultraviolet becomes necessary as the red- 
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Figure 4. Zoom into the merging zone for MILES-based (green) 
and UVBLUE-based (black) models, for various ages and metal- 
licities. 



shift gets progressively higher. Of the four sets of 
models, only the Pickles-based has a UV-extension. In 
order to improve on this situation, alternative ver- 
sions of models have been constructed, where we have 
performed a c o ncaten ation with the SSP models of 
iMaraston et al.l (|2009al ) ■ which are based on the fully 
theoretical UVBLUE/Kur ucz library of high-resolutio n 
(R= 10000) stellar spectra i|Rodriguez-Merino et ail 120051 ). 
and have a wavelength coverage of 1000-4700 A. After down- 
grading the UVBLUE-based models to the appropriate res- 
olution of each respective set of models, and normalising 
them onto the M05 models using the average flux in a 100 
A wide passband centered at 3050 A, the merging takes place 
around 3750 A. For the ELODIE models the merging is per- 
formed between 3960 and 3990 AQ 

Figure [4] zooms into the merging zone between 
UVBLUE-based and empirically-based MILES models for 
various ages and metallicities. The overlapping region looks 
excellent at young ages, where also the population flux in 
the UV is higher. At intermediate-ages the matching is still 
very good, though one notices some excess flux in the the- 
oretical library especially at super-solar metallicity. At old 
ages the excess flux is evident at all metallicities. It is hard 
to say what the source of this noise is, but the flux level here 
is very low and should be irrelevant in spectral analysis. Fi- 
nally, Figure [5] shows another example of merged models, us- 
ing the 400 Myr STELIB-based model (blue) and the M05. 
The merging takes place here at 3752.5 A. 



1 Attempting to keep the empirical approach intact, we also tried 
to take complementary stars from the Pickles library. We found, 
however, that the number of lower resolution Pickles stars needed 
to appropriately sample the lower MS, was enough to affect the 
overall resolution of the STELIB-based SEDs. 



2 The exact wavelength changes slightly for aesthetical pur- 
pouses. 
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Figure 5. Comparison between the merged UBVLUE - STELIB 
400 Myr SSP (blue) and the equivalent M05 model (black), and 
zoom into the merging region (inset). 



4.6 High-resolution theoretical SSPs based on the 
MARCS library 

High-resolution models that reach further into the near- 
infrared than what can be accomplished with empirical spec- 
tra are equally useful. We therefore calculated a grid of SSP 
models adopting the M ARCS library of high-r esolution, the- 
oretical stellar spectra l|Gustafsson et al.l2008h . which covers 
an impressive spectral range from 1300 A to 20 /xm with a 
resolution of R=20,000. The library is somewhat limited in 
terms of stellar temperatures - the highest temperature is 
8000 K at any metallicity - which translates into a minimum 
age of 0.8 Gyr at 2 Zq , and around 6 Gyr at the lowest metal- 
licity. In the latter case a problem arises also for the oldest 
populations, where the shedding of mass towards the end 
of the RGB phase causes deeper, hotter layers in the stars 
to be exposed. If the mass-loss is severe enough, horizontal 
branch stars hotter than 8000 K will be produced. For a de- 
tailed account of the mass-loss prescription employed here, 
the reader is referred to M05. 

Furthermore, we calculated a version of the MILES- 
based SSP models in which we extended them into the 
UV as explained in Section 4.5 and into the near-IR using 
the MARCS-based models illustrated in this Section. These 
merged models are probably the models with the largest 
wavelength coverage at high resolution that are presently 
available. 




log (Age/yr) 

Figure 6. Schematic view of the available Pickles-based models, 
in terms of age and metallicity, Z. 



5 THE SSP MODELS 

5.1 Age and metallicity coverage as a function of 
the spectral library 

The simple stellar population (SSP) models that can be cal- 
culated for each of the four different empirical libraries, in 
terms of age and metallicity, are shown schematically in Fig- 
ures 1 61 9 1 It is evident that empirical spectra alone cannot ac- 
count for all desired combinations of these parameters. For 
example, no library contains the hot stars of low metallic- 
ity that are required to build metal-poor SSP models with 
young ages. Naturally, this is because star formation in the 
Milky Way takes place in the metal-rich disc. 

As mentioned earlier, the Pickles-based SSP models 
(Figure [6} are only available at solar metallicity. The pres- 
ence of O-type stars in the library allows for very young 
ages to be modelled (down to around 2 Myr, where the 
turn-off stars have a temperature of roughly 40000 K), but 
these spectra are also significantly attenuated by dust ab- 
sorption in the UV. Hence the Pickles-based models with the 
youngest ages (below ~100 Myr) should be used somewhat 
more cautiously. 

The MILES-based SSP models display the most com- 
prehensive sampling of the model grid of all the four empiri- 
cal libraries (see FigQ . Reflecting the chemical evolution of 
the Milky Way, no particularly young ages can be modelled 
at the lowest metallicity. In general, the mapping of stellar 
atmospheric parameter at this metallicity is rather coarse, 
hence these models should be treated with some reserva- 
tion. As previously mentioned, the [Fe/H]= —2.0 models, 
together with the [Fe/H]= —1.3 models, are expected to be 
enhanced in [Mg/Fe] by around +0.4 dex. 

The amount of spectra in the STELIB library with twice 
solar, solar, and half solar metallicities is adequate to pro- 
duce models down to relatively young ages. A caveat is that 
some of the stars of non-solar metallicity in the STELIB 
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Figure 7. Same as Fig[6]for the MILES library. 



Figure 9. Same as Fig|6]for the ELODIE library. 
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Figure 8. Same as Fig[6]for the STELIB library. Models in the 
striped area have limited wavelength coverage. 



library do not have complete spectroscopic observations 
over the entire quoted wavelength rang e of 3200 - 9300 A. 
The correspondent models will cover a narrower wavelength 
range of 3200 - 7900 A (see Figure [8]). 

The ELODIE library contains a number of very hot 
O-type stars, which in principle would allow us to produce 
models as young as 1 Myr. The spectra of these recently 
born stars are, regrettably but not surprisingly, severely red- 
dened by close-in dust, and/or displaying various amounts 
of emission lines, and many of them had to be discarded 
in the population synthesis, thereby setting the age limit at 



solar metallicity to ~3 Myr. At half-solar metallicity it is 
instead a lack of hotter supergiants (> 6000 K) that sets a 
limit to the age of the youngest model. 

It should be noted, though, that despite the removal of 
the worst spectra, the youngest models are not completely 
unaffected by the dust reddening of their turnoff stars. At 
the lowest metallicity, 1/100 Zq, the sampling of (T e //, 
log(g))-space in the ELODIE library is fairly crude. A few 
ages can be modelled, but a complete evolutionary coverage 
is lacking, and they should thus be treated with some cau- 
tion. There are not enough stars around [Fe/H]= —1.30 to 
warrant the calculation of models of any age at that metal- 
licity. 



5.2 The SSP spectra 

5.2.1 MILES models at solar metallicity 

Figure [10] shows the age effects at fixed solar metallic- 
ity, comparing our M ILES-based models to the models by 
IVazdekis et~al1 (|201Ch (that are based on MILES). The 0.1 
Gyr model compares well. The higher flux longward 5000 
A in the models at 1 Gyr is due to the TP-AGB phase. At 
old ages, the slightly higher flux longward 6000 A in the 
Vazdekis models is due to the redder RGBs of the Padova 
tracks that are used in these models (see Figure 7 in M05 
and discussion). 

We also calculated broad-band colours, in order to spot 
differences due to the different libraries. 



5.2.2 Young models as a function of spectral library 

Figure [11] displays the integrated colours of the youngest 
models as a function of the various spectral libraries. Colours 
plotted are those that could be computed given the intrinsic 
wavelength extension of the library through the near-IR, e.g. 
MILES allows the computation of the sole B — V, whereas 
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Figure 10. MILES-based models (green) for various ages and a solar metallicity, compared to models by Vazdekis et al. (2010, magenta). 
The latter models are not provided for 10 Myr of age. The higher flux longward of 5000 A in our 1 Gyr model is due to the TP-AGB. 




log Age (yr) log Age (yr) 

Figure 11. Integrated colours of young SSPs as a function of the various empirical spectral libraries. Colours plotted are those that 
could be computed given the wavelength extension of the library, see the text. The Pickles-based models using the M05-UVBLUE run 
identical and have been omitted. 
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Figure 14. Effect of changing the temperature scale on stel- 
lar population models, for the illustrative example of a 12 Gyr 
solar metallicity SSP. STELIB stars have been tied onto the 
MILES temperature scale by using the fit of Figure 3. The original 
STELIB-based model (blue line) turns into one with lower flux 
in the near-IR (black line) when the MILES temperature scale 
is adopted. For comparison, the original MILES-based model is 
also plotted (green). The near-IR slope of the models depends 
substantially on the adopted temperature scale of RGB bump 
stars. 

STELIB allows also the computation of V — R and V — I. 
U — V and U — B are calculated through the extension with 
the M05-UVBLUE models. The Pickles-based models based 
on their intrinsic UV from Pickles (which in turn is based 
on empirical spectra from the IUE database, cf. Section 3.1) 
and the version extended with the UVBLUE library run very 
similar, hence we show only the first ones. 

Overall, one sees an excellent agreement in the colours 
of very young stellar populations as a function of the spectral 
library for most colours. Some discrepancy is visible in the 
optical-to-near-IR colours of model ages where red super- 
giants exist (t ~ 10 Myr), which are however well inside the 
typical scatter associated to colour indices at these young 
ages, which in turn is due to the short duration of the red 
supergiant phase (see Marigo et al. 2008 for a recent plot; 
see also Langon et al. 2008 for a general discussion on red 
supergiants). 

Note that the good matching for the continuum blue- 
ward of 4000 A supports our method of merging the 
empirically-based models with the theoretical UVBLUE- 
based models. 

5.2.3 Old models as a function of spectral library 

Figure [12] displays the overall spectral behaviour of the vari- 
ous models for an old, solar metallicity population. For clar- 
ity models have been smoothed to the quoted Pickles reso- 
lution. M05 is shown as a black line. 



In general, significant differences are found between 
models based on different empirical spectral libraries. 

Shortward A ~ 5000 A (top panel) the models run very 
similar and close to the M05 model, the latter displaying a 
small excess flux around 4600 A. In the range 5200-6400 A 
the various models display a very clear behaviour, i.e. all 
empirical models display a flux deficit with respect to the 
M05. We focus on this effect in the next subsection. The 
only exception to this trend are the ELODIE-based models, 
which have an overall higher flux, lying in between the other 
models and the M05, and reconciling onto the M05 at A > 
6300. The reasons for such trend remain unclear, but it could 
possibly be related to the flux calibration, which at very high 
resolution is quite difficult. 

In addition, as noticed in iMaraston et alj (|2009al ) and 
discussed in Section [31 the MILES spectra lack flux long- 
ward of A ~ 6400 A which is reflected in the SSP models 
(Figure HU). 

We find that a substantial part of this effect is due to 
the assumed temperature scale for RGB-bump stars with 
temperatures between 4000 and 5500 K. In Section 3.6 we 
discuss the difference in temperature scales and we made 
an experiment in which we tied the STELIB library onto 
the MILES temperature scale. Here we show the resulting 
stellar population models in Figure ITU The blue line shows 
the original STELIB-based SSP, the black line the one we 
obtain using the MILES temperature scale, and the green 
line shows the MILES-based SSP. It is striking how much the 
near-IR flux has dropped by applying the simple re-scaling. 

We could figure out that this effect is dominated by 
RGB-bump stars. 

This highlights the importance of understanding the ef- 
fects of the assumptions o n adopted temperature scales on 
stellar population mo dels dPercival fc Salarijl2009l) . For ex- 
ample, as discussed in IMaraston et al.l (|2009a), the different 
slopes of the models affect optical-to-near-IR colours and 
the interpretation of galaxy data. 

We notice that our models based on the Pickles and 
STELIB libraries, but also those based on the theoretical 
library MARCS, in which no temperature scale calibration 
needed to be input, agree well long-ward 6000 A (cfr. Figure 
15). Given the effect a different temperature scale for RGB 
stars has on the spectral shape of the integrated models, in 
Appendix [X] we further expand on this point and provide 
a version of our MILES-based models in which we revise 
the near-IR slope to better match the one of the other li- 
braries. Though we cannot be certain that the near-IR slope 
of MILES is not correct only because it does not compare 
well with those from other libraries, another fact in support 
of a revised slope is the observed-frame r — i colours of mas- 
sive luminous and red galaxies in the 2SLAQ survey, which 
were found to be systematically redder than colours ob- 
tained with a variety of existing models (Wake et al. 2006). 
In Maraston et al. (2009) we improved upon that issue us- 
ing Pickles-based models, and we noticed that MILES-based 
models - with their intrinsic near-IR slope - did not perform 
as well (Figure 3,4 of Maraston et al. 2009). Though not 
being a clear-cut argument, it is an interesting hint towards 
a specific spectral shape. Obviously, observations of simpler 
systems such as star clusters would be very beneficial to 
constrain the effect. 

Finally, at A ~ 6500 A (middle and bottom panels) 
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Figure 12. Comparison of the four sets of high-resolution empirical models and the original M05 models, normalised 
to 5050 A. SEDs refer to 12 Gyr, solar metallicity SSPs with RHB and a Salpeter IMF. Models have been smoothed 
to the quoted Pickles resolution. 



Pickles-based and STELIB-based models, and M05 run quite 
similar one to each other, except at 7500 A 8300 A and 9000 
A where the Pickles-based models have consistently lower 
flux than the M05 models, whereas the STELIB-based mod- 
els lie closer to M05. It is presently quite hard to calibrate 
these effects as one would need excellent quality near-IR 
spectra of old, solar metallicity GCs such as those in the 
Baade window, that are on the other hand strongly affected 
by reddening. While these will be available in the near fu- 
ture, we shall use galaxies to try discriminating among these 
effects in a future paper. 

It should be stressed that in Figure [12] we have nor- 
malised all models at a certain wavelength (5050 A), which 
allows us to emphasise the spectral shape around it. This 
is obviously just a visualisation, the quantitative effect is 
derived by calculating the flux ratios of the models, i.e. by 
means of broadband colours. These are shown in Figure [131 
which is analogous to Figure [11] for older ages. 

The colours of old models are contributed to by several 
stellar evolutionary phases (see e.g. M05) and the compari- 



son of Figure [13] tells us that - at the same energetics - dif- 
ferences in the individ ual stellar spectra matte r. This point 
was amply discussed in lMaraston et all l|2009al ) where a first 
version of these models was published, and in M05 regarding 
MW GCs. 

In particular, we confirm the results of iMaraston et all 

(2009a|) for the g' — r' colour, here for the B — V colour, which 
when using empirical libraries - all of them - is bluer than 
the one relative to the M05 models, which are based on the 
BaSeL assembly of the Kuru cz model atmospheres . In par- 
ticular, notice that both the IVazdekis et~ai1 (|201Ch models 
and the M05-MARCS feature a bluer B — V, in close agree- 
ment with the empirical models. This point will be amply 
discussed in the next subsection. 



5.3 The spectral shape around the V-band 

In lMaraston et ail (|2009al ) we found that SSP models based 
on empirical libraries (from Pickles, MILES, STELIB) dis- 
played a lower flux around the l/-band with respect to 
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Figure 13. Integrated colours of old population models as a function of the various spectral libraries. The red dashed line refers to 
Pickles-based models in which their UV spectrum has been substituted with the theoretical UVBLUE spectrum (see Section 4.5). The 
models by Vazdekis et al. (2010) are also included. 



equivalent models based on the BaSeL theoretical library, 
an effect which helped to a better matching of the SDSS 
g',r',i' colours of Luminous Red Galaxies around redshift 
0.4. There, we also showed (Figure 4) that the same effect 
was found in the individual stars from the libraries, espe- 
cially in giants stars, which secured us against possible arti- 
ficial effects introduced by the synthesis itself. The models 
used in that paper were a sub-sample of the extensive grid 
presented here, where we confirm our earlier finding. 

In this paper we consider two further models that 
stre ngthen our conclusions. F igure ll5l is equivalent to Figure 
3 in iMaraston et al.l l|2009al ) , and a blow-up of the middle 
panel of Figure [T2l with two important addition, namely the 
IVazdekis et alj (|201Cl ) models based on MILES (magenta) 
and our new models based on the theoretical stellar library 
MARCS (Section 4.6). The figure clea r ly sh ows that the 
completely independent IVazdekis et al.l (|2010T ) models dis- 
play exactly the same effect as ours. Furthermore, the M05- 
MARCS models behave as the empirically-based models. 
Note that the latter models do not require any adjustment 

3 With the exception of the ELODIE-based models as noticed 
earlier 



or scaling of the MARCS spectra to be computed, which 
reinforces the conclusion that the procedure we adopt to in- 
sert the empirical libraries in the synthesis is appropriate. 
Finally, as we already pointed out, the slope long- ward 6000 
A agrees very well between STELIB-based, Pickles-based, 
and MARCS-based models, while MILES-based models dis- 
play lower flux, which is mostly due to the temperature cal- 
ibration of RGB stars in this library. 

For completeness, we also show in Figurc[T|J]the individ- 
ual stellar spectra from the various librarie s, including the 
MARC S library, in a remake of Figure 4 of IMaraston et all 
( 20093). The MARCS synthetic stellar spectra behave as the 
empirical libraries at given stellar parameters. Finally, Fig- 
ure [T7] shows that the B — V colour of the individual stars 
from the empirical or the MARCS library is bluer than those 
of the Lejeune stars. 

Note that this result is obtained by trusting the stellar 
parameters provided by the libraries. The agreement with 
the MARCS-based modes suggests that this is the right pro- 
cedure to adopt. 

Figure [TH] shows the models of IConrov fe Gunnl (|2010T ) . 
These models behave as ours and the Vazdekis et al. models. 
This is at odds with the conclusions of Conroy and Gunn 
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Figure 15. Remake of Figure 3 of Maraston ct al. (2009b). The plot compares the spectral region around the V-band for different 
models, all referred to a 10 Gyr SSP with solar metallicity and a Salpeter IMF. New here is the inclusion of the models by Vazdekis 
et al. (2010) based on MILES (magenta line) and models based on the theoretical library of stellar spectra MARCS (orange line). All 
high-res. models have been smoothed to the Pickles's resolution. 
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Figure 17. B — V colours of individual stellar spectra from vari- 
ous libraries. The aim of this plot is to compare the optical colour 
between empirical libraries - using MILES as representative of 
them - and two theoretical libraries, namely the Lejeune et al. 
compilation of the Kurucz library and the MARCS library. 



who cl aim to not recover the findings of Ma raston et al.l 
l|2009al ), i.e., they were unable to find any difference between 
the models based on the theoretical BaSeL library and the 
empirical MILES library, in particular, not on the B — V 
colour. Figure [19] shows the B — V colour of their models, 
based on MILES and on BaSeL (red dotted and solid). Op- 
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Figure 18. As Figure fl5l for the lConrov fc Gunnl lj2010l ) models 
(for a Chabrier IMF). 



posite to their statement, the colours are different, and in 
particular those based on the MILES library are bluer, as in 
our models (black lines in the same plot). 

Conr oy fc Gunn (2010) spec ulate that the effect de- 
scribed in lMaraston et all (1200 9b) is due to the "particular 
way in which the empirical stellar spectra were inserted in 
the synthesis by Maraston et al. (2009b)" but as we clearly 
show in this paper, their models - in spectra and in colours 
- behave the same way as ours. 

Conroy & Gunn (2010) also cite Bruzual & Chariot 
(2003) for having reached their same conclusions on the 
B — V colour. Instead, we show in Figure [T9"l that the Bruzual 
& Chariot (2003) models - at odds with the cite by Conroy 
& Gunn (2010) - also display an effect due to the stellar li- 
brary, which has however the opposite sign, i.e. the STELIB- 
based colours are redder than those based on BaSeL. Here we 
note that these authors adopt a different procedure to insert 
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Figure 16 . Individual stel l ar spec tra from various libraries for several sets of stellar parameters (labelled in the panels). Remake of 
Figure 4 of Maraston et al.l ll2009al ) including now the MARCS library. 
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Figure 19. Integrated B — V colour as a function of the spec- 
tral library, for different models, by Conroy & Gunn (2010, red); 
Bruzual & Chariot (2003, blue); M05 and this work (black). Solid 
lines refer to models based on empircal libraries, dotted to those 
based on the BaSeL library. The Conroy &: Gunn models refer to 
a Chabrier IMF, but this implies a minimum effect. 



gin of the effect, but is beyond the scope of this paper to 
investigate this issue in more depth. 

5.4 The near-IR spectrum of old populations 

The Pickles-based models extend to the Jf-band, hence en- 
abling us to check the near-IR wavelength regime of the 
spectra as a function of stellar library. We shall compare 
these with two variants of the M05 models based on theo- 
retical libraries, the published M05 that was based on the 
BaSeL-Kurucz library, and the M05-MARCS that we pub- 
lish herfl Figure [20] shows the SEDs of 10 Gyr old pop- 
ulations according to the various models (we checked that 
the main features are invariant with the age). Though the 
overall shape is consistent, there are two regions where the 
models deviate from each other, namely around 1 micron, 
in correspondence to the J-band and close to 1.5 micron, in 
correspondence to the //-band. The behaviour is not sim- 
ple, the models cross, with the standard M05 having more 
flux in the J-band and less around H, and viceversa for 



the empirical spectra in the population synthesis, which is 
matching the empirical spectra to the theoretical ones (from 
BaSeL) in the B — V colour, i.e. without using the library 
stellar parameters. This different procedure may be the ori- 



4 Notice that both these models based on theoretical spectra ex- 
tend to much larger wavelengths, but here we are interested in the 
main emission from stars and in comparing with the Pickles-based 
models so we plot just up to the i\-band 
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Figure 20. The near-IR spectrum of various SSP models, the one 
based on the empirical Pickles library (red), and two based on 
theoretical atmospheres, the standard M05 based on the BaSeL- 
Kurucz library (black) and the one presented here based on the 
MARCS theoretical spectra (orange). Models are 10 Gyr old and 
have solar metallicity. 



M05-MARCS and Pickles-based. These trends explain the 
difference in the corresponding broad-band colours. From 
the observational side the near-IR spectrum of galaxies and 
star clusters is becoming a new promising area, of research 



(e.g.jRiffel et alj|2007t 120081; iLancon et al. pool; ISilva et al 



20081 ; ICesetti et al.ll2009l ; lLvubenova et aLlkoialRiffel et a! 
201ll ; lMiner et al.ll201ll ) thanks to new instrumentation such 



as XShooter and others, and in preparation for the James 
Webb Space telescope. We shall soon be able to check these 
spectral details. 



5.5 Lick indices 

A proper calculation and analysis of Lick absorption in- 
dices in integrated stellar populations require the account 
of element abundance ratios, because these lines are sensi- 
tive to their pattern, as amply discussed in the literature 
ilWorthev et all Il992l; bayies et al.lll993l; I Surma fc Benderl 



1995 
2000 



Tripicco fc Belli Il995l; iGreggid 1 19971: 



Maraston ct al. 



2003 



g| Thomas et all 



Schiavonl 120071 ; iThomas et al-lHoill , .iust to 



Trager et al. 
2004: 



2003 



name a few 

among the most relevant works on the topic, see the lat- 
ter paper for full referencing and discussion). 

In particular, Figure 1 of lMaraston et all <|2003l ) clearly 
shows that the Mgb and Fe Lick lines of metal-rich galactic 
Bulge globular clusters with measured [a/Fe] ~ 0.3 are not 
matched by stellar population models adopting the stellar 
index line-strengths of Milky Way solar-neighborhood stars, 
because the latter have a total solar metallicity and solar- 
scaled element ratios while the former have the same total 
metallicity, but [a/Fe] enhanced ratios. Star clusters in the 
Milky Way with abundance ratios known from individual 
stellar spectroscopy are not as ambiguous as distant galaxies 
and were used by Maraston et al. (2003) to definitively unveil 
the nature of the persistent discrepancy between Lick mod- 
els and galaxy data. The models by Thomas et al. (2003), 
that account explicitly for various elemental ratios (see the 
same Figure), match well the Bulge globular cluster indices 
for [a/Fe] ~ 0.3, consistent with the measurements from 



integrated-light spectroscopy. They also match the lines of 
massive galaxies. 

The MILES library is not much different from the early 
Lick library from the point of view of abundance-ratios, in 
that it contains the same Milky Way-type of stars and in 
particular no Bulge stars, that would provide the metal-rich 
and [a/_Fe]-enhanced population that is missing in the solar 
vicinity. 

Hence, the indices calculated on the model SEDs we 
present in this paper are expected not to match the lines 
of metal-rich and [cv/_Fe]-enhanced populations, that can be 
found i n external massive ga laxies and star clusters. For this 
reason. pThomas et al.l (|201ll ) calculated element-abundance- 
ratio-sensitive models based on MILES stars but following 
the same philosophy of Thomas et al. (2003), i.e. by theo- 
retically manipulating the index response functions. 

Nonetheless it is interesting to see how the Lick indices 
calculated on the integrated model SEDs of this paper com- 
pare with those from the TMJ models. Figure [2T1 shows all 
25 Lick indices in order of increasing wavelength from top 
to bottom, with the exception of the Balmer lines that are 
listed besides one each other at the top. Lick indices di- 
rectly calculated on the model SEDs presented in this pa- 
per are shown as black lines, TMJ models for two differ- 
ent abundance-ratio patterns, namely solar-scaled (red) and 
[a/Fe] = 0.3 (blue), are also shown, with solid and dashed 
line-style referring to solar and sub-solar (~ —1.35) total 
metallicity [Z/H]. Also shown as green lines are models ob- 
tained by integrating the MILES-based FF of Johansson 
et al. (2010). This type of models help us understanding 
whether discrepancies between the SED-based and the TMJ 
models are due to the SED vs FF approach. 

In general, indices from the integrated SEDs (black 
lines) should be in agreement with solar-scaled TMJ models 
at solar metallicity (solid red lines) and with [a/Fe] — 0.3- 
TMJ models at low metallicity (dotted blue lines), because 
of the intrinsic chemical bias Milky Way stars carry into the 
calculations (see discussion above). This is indeed the case 
for the vast majority of indices. Only a few cases do not 
follow this trend, on which we comment below. 

As mentioned, the comparison between green lines and 
black lines helps us spotting discrepancies due to the differ- 
ent approach, namely integration on SEDs or FF. Overall, 
there is an amazing agreement between our calculations fol- 
lowing the two approaches for most indices at both metallic- 
ities. Neglecting minor differences, such as e.g. the Ca4227, 
one index sticks out for having a discrepant value between 
direct integration (black) and FF approach (green), namely 
Fe4531 at sub-solar metallicity. We verified that this is due 
to a slight offset of the correspondent FF for giant stars. The 
Fe4531 indices of these stars are particularly sensitive to 
gravity effects (Johansson et al. 2010). Thus the modeling 
of this index is complicated, from both the SED and FF 
point of view, especially at metallicities below —1.0 where 
the number of MILES stars is fairly low. However, the chem- 
ical enrichment modelling by TMJ is able to rectify the off- 
set introduced by the FF, such that the end-point TMJ a- 
enhanced model (blue line) perfectly agrees with the values 
integrated onto the SEDs (black lines), and the TMJ index 
model is perfectly calibrated against GCs (see TMJ, Fig- 
ure 3). 

At sub-solar metallicity some indices do not follow the 
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Figure 21. Comparison between different model Lick indices (at Lick resolution), namely: directly calculated on the MILES- based in- 
tegrated SEDs presented in this paper (black lines); calculated through MILES-based FF (from Johansson et al. 2010, green line); from 
models with variable abundance ratios from TMJ for solar-scaled [a/Fe] = 0.0 (red) and [a/Fe] = 0.3 (blue). Solid and dashed line-styles 
refer to solar and sub-solar ~ — 1.3 metallicity. Vazdekis et al. (2010) models also directly calculated on SEDs and for the same metallicities 
are shown as magenta lines. 



expected trend of matching the TMJ [a/_Fe]-enhanced (dot- 
ted blue lines). These are: Hf, which is weaker in the SED- 
based approach with respect to the TMJ model, and the 
G-band G4300 - and to a lesser extent Fe4383 and Fe5406 - 
which are stronger. In all these cases one sees that the use 
of FF in TMJ vs the SED-integration approach (i.e. green 
vs black) is not responsible for the effect (with the possible 
exception of H'y at low ages, where the warmest MS dwarfs 



have a higher index value than the one approximated by the 
FFs). The major part of the effect comes from the correction 
of the element ratio bias in the TMJ approach. As discussed 
in TMB, these models assume that the stellar 'metallici- 
ties' [Fe/H] in the input library reflect true iron abundance 
rather than total metallicity. As a consequence, the metal- 
licity scale gets slightly shifted through the bias correction, 
which leads to weaker metal and stronger Balmer indices at 
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low metallicities. This effect is indeed seen here with G4300 
and Hy being most affected such that the index G4300 is 
weaker and H-y stronger in the TMJ model. It should be 
noted that these particular indices follow a very good cali- 
bration with GC values as shown in TMJ. 

At solar metallicity, there is none of such noticeable 
cases, and the SED-based calculations (solid black) and the 
TMJ models with solar-scaled abundance ratios (solid red 
lines) agree well. Some slight differences can be seen for 
Ca4227, C 2 4668 and CN1, which ought to be ascribed to 
the FF being slightly lower than the SED calculation. 

We note that the indices by Vazdekis et al. (2010) (ma- 
genta lines in Figure I2ip - which are obtained through di- 
rect integration on their MILES-based model SEDs - display 
similar discrepancies for the above mentioned indices, but 
these are typically larger and show a different age depen- 
dence, see the G4300 band, 6*24668 or the Fe5015, among 
the others. Note also that some indices are well offset, such 
as the Ca4227 and the TiO indices. The exact reason for 
these difference is hard to assess, but it should be kept in 
mind that the Vazdekis et al. (2010) models are based on 
different evolutionary tracks (from Padova) which have red- 
der RGBs than the tracks used in our models and a different 
overshooting (see M05). 

Aside from the above discussion, note that the TiO in- 
dices are stronger at ages around 1 Gyr due to the TP-AGB 
contribution in the Maraston models and the use of the em- 
pirical TP-AGB spectra (see Section 3. 7)0. Also the Mgj 
index, which is known to be sensitive to Carbon (Thomas, 
Johansson, Maraston 2011; see also Tripicco & Bell 1995), 
is stronger at the same ages in the SED-based models. 



6 TESTING THE MODELS WITH GLOBULAR 
CLUSTERS 

In this section we provide a few examples on how the new 
models perform in fitting real objects. Following our pre- 
vious works, we use globular clusters as they are the clos- 
est approximation to simple stellar populations, i.e. age and 
metallicity for each cluster can be determined from resolved 
studies independently of stellar population models. Suitable 
for this test are Milky Way and Magellanic Clouds objects. 
As we are mostly interested to test the effect of spectral li- 
braries, which - as they do not contain any TP-AGB star 
- is mostly concentrated on old populations, we shall not 
repeat all comparison plots published in M05 using Magel- 
lanic Clouds clusters. In the following we shall use only MW 
GCs data. As is well known, MW field and globular cluster 
stars feature a well defined pattern of chemical abundance 
ratios, in particular the [a/Fe] ratio - at given total metal- 
licity. This introduces a bias in models based on empirical 
MW stars (Maraston et al. 2003; Thomas, Maraston & Ben- 
der 2003) that we shall take qualitatively into account in our 
comparisons. Note that empirical libraries contain - besides 
field stars - also globular cluster stars, in particular bright 
red giant branch stars, but these show a similar element 



5 The MILES library docs not contain TP-AGB star spectra 
hence these are not included in the TMJ fitting-function-based 
models 



abundance pattern as field stars (|Pritzl et al.ll2005h hence 
their inclusion does not affect our conclusions. 

In the following we show tests involving both GC spec- 
tra as well as GCs broad-band colours. 



6.1 Globular cluster spectra 



We use the samples of iPuzia et al. (2002) and 
Schiavon et all l|2005f ). with main focus on the latter, 
which consists of 40 objects that represent the full range of 
MW GC parameters at an instrumental resolut ion of 3.1 A 
FWHM, and S/N ratios between 100 and 400. IPuzia et all 
(2002)'s spectra have a coarser resolution (6.7 A FWHM) 
because the aim was to evaluate Lick indices, and have 
concentrated on the most metal-rich Bulge clusters, but 11 
out of their total number of 12 objects are in common with 
the lSchiavon et all (|2005h compilation, thus enabling a com- 
parison also between different samples (see Section 6.1.3). 
Both sets of spectra have been flux calibrated, but neither 
have been corrected for reddening. 

GC ages and metallicities from our new SED models 
have been derived through a standard SED-fitting proce- 
dure, where we compute the x 2 between observed and model 
SEDs over the wavelength range 3600-6200 A, pixel by pixel 
(2600 pixels), using the S/N spectra for assigning errors to 
each pixel[f] The observed spectra are de-reddened prior to 
the fitting, using a simplified version of the wavelength- 
depen dent MW interstellar extinction law of iFitzpatrickl 
(1 1999T) with R=3. 1 and the E(B — V) parameter compiled 



Schiavon et ail (2005) This parameter depends strongly 



on the line-of-sight; however, since the purpose of this exer- 
cise is to show that, despite a rather simple approach, very 
reasonable estimates of GC parameters can be obtained in 
all but a few cases, such approximations may be forgiven. 

For testing we use both the MILES-based and STELIB- 
based models, which have a resolution comparable to the GC 
spectra. Only the MILES-based models could be used for all 
GCs, as they exhibit a wide enough metallicity coverage (cf. 
Figure 7). Model grids have been interpolated in age and 
metallicity in order to avoid clustering of solutions around 
the grid points. The Horizontal Branch (HB) morphology 
has been taken into account, in the sense that model SEDs 
have been pre-selected according to ifQ- SED fits have been 
performed using both Salpeter and Kroupa (2001) IMF- 
mo dels. Results wi ll be c omp ared to CMD-derive d ages , 



by |Pe Angeli et all (|2005h and iMarfn-Franch et al 



2009) 



1996) 



and metallicities from the latest version of the lHarris! 
catalogue of Milky Way Globular Clusters and from Marin- 
Franch et al. (2009). 

Figure |2"2"1 shows examples of SED fits of GCs, with the 
best-fitting MILES-based models shown in red. The six clus- 
ters represent a fair sampling of the MW GC parameter 
space, and are also those with the least amount of redden- 
ing (E(B-V)^ 0.06). In general, the fits are excellent and 



6 The Schiavon et al. spectra are corrupted around 4546 A and 
5015 A and this range has been masked out. 

7 The models of this paper as the M05 models are provided for 
two morphologies of the Horizontal Branch, namely intermediate- 
red and blue, following a standard recipe for mass-loss that was 
calibrated on Milky Way GCs. See M05 for details. 
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the parameters of the best fits are in Rood agreement with 
the age values d erived i n De Angeli et al.l l|2005l ) and the 
metallicities from lHarrisl lfl996fl. 

As pointed out by Vazdekis et al. (2010), the discrep- 
ancy around 3850 A visible for NGC 104 (47 Tuc) is due to 
the strong CN features of t his cluster compared to the stars 
composing the libraries. In lThomas et al.l l|201l[ ) it is shown 
that the corresponding CN Lick indices of such GCs can be 
matched by enhancing Nitrogen. 

We also spot a recurring discrepancy between 4000 and 
4100 A, where the model flux is slightly weaker compared to 
the observed spectra. Finally, the models display somewhat 
lower flux levels blueward of ~ 3750 A, but here the S/N has 
decreased roughly a factor 3 compared to the values around 
5000 A. A similar comparison is provided by Vazdekis et al. 
(2010) for the clusters NGC 104, 5904 and 7089 (their Fig- 
ure 20), with derived cluster parameters in good agreement 
with ours (also for NGC 7089, not shown here). 

So far we have focussed on the comparison of the over- 
all SED. We do not attempt a detailed comparison of indi- 
vidual absorption features, because the latter requires the 
account of element abundance ratios in the models, which 
is done in a series of companion paperfl Thomas, Maras- 
ton & Johansson (2011, hereafter TMJ11) present the flux- 
calibrated version of the widely-used Thomas, Maraston & 
Bender (2003) model Lick indices with variable abundance 
ratios. These have been upgraded using empirical calibra- 
tions of the index stre ngth as a function of s tellar parame- 
ters (fitting fun ctions, Johansson et al.ll2010T l based on the 
MILES library. iThomas et all (|201lh "address in detail the 
comparison of these latest model Lick indices with Milky 
Way globular cluster data. 

What is interesting to evaluate here is - however - which 
ages and metallicities we get for MW GCs had we used - 
instead of the full SED spectral fitting - the Lick indices 
calculated on the same SED models (see Section |5.5[1 . i.e. 
without an explicit account of element abundance ratios. 
These are shown as filled symbols in the central panels of 

Figures [23 [sa rjzirji 

Ages are strongly underestimated with respect to the 
CMD ages (y-scale), which is due to the metallicity being 
overestimated (Figure [2U central panel). The situation im- 
proves when considering a small subsample of indices (open 
symbols), the combination of which is little sensitive to ele- 
ment abundance ratios (namely, Hg A , Mgb, Fe4383, Fe5335, 



8 In considering these comparisons it is important to remember 
that we show the peak values in the marginalised \ 2 distributions. 
Using a 68 % confidence interval we obtain very large age and 
metallicity ranges - a perfect showcase for the age-metallicity 
degeneracy. On the other hand, in most cases the peak values are 
in excellent agreement with the CMD determinations. Two such 
examples are given in Figures 1231 and 1241 where the normalised 
age- and metallicity probability distributions are highlighted for 
one of the best (NGC 1851) and one of the worst (NGC 6544) fits. 
For NGC 1851 we obtain 9.0 Gyr and [Fe/H]=-1.0 vs CMD values 
of 9.0 Gyr and [Fe/H]=-1.2, with a confidence interval of 0.7-15 
Gyr, and -1.7 to +0.3. For the highly-reddened cluster NGC 6544 
literature values are 9.4 Gyr and [Fe/H]=-1.56, and we get 8.0 
Gyr and [Fe/H]=+0.3 with a 68% confidence interval of 1.5-15 
Gyr, and -1.0 to +0.3. 

9 Spectral features at other wavelengths for the present models 
will be addressed in future works. 
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Fi gure 23. Marginalized age- and metallicity \ distributions for 
the iSchiavon et al.l J2005h globular cluster NGC 1851 from fitting 



with MILES-based models, normalized at the peak values. This 
is one of the better fits. 



Fe5406). This suitable index combination was defined by 
TJM11 as a base index set. 

Left-hand panels in the same figures shows the peak val- 
ues of the marginalised age and metallicity \ 2 distributions 
derived through full SED fit (cf. Figure 22). Right-hand pan- 
els show the values of age and total metallicity derived by 
IThomas et all l|201ll) based on their model Lick absorption 
line indices. GCs ages, metallicities and [a/Fe]-enhancement 
are derived through a x 2 -fitting using the best-calibrated in- 
dices (see paper for details). 

Ages derived with both a full SED spectral fitting 
or using selected spectral lines with detailed modelling of 
abundance ratios are in much better agreement with the 
CMD-derived ones compared to those obtained using in- 
dices without a modelling of abundance ratios, even if the 
CMD-derived ages span a narrower range than those from 
integrated-light methods^ This gives us confidence in the 
intrinsic age-scale of the models. 

Relevant to this test is to note that CMD-derived ages 
carry their own problems as there are zero-point issues de- 
pending on the adopted tracks and on whether element ra- 
tios are taken into account. Figure [27] shows an age compar- 
ison identical to Figure 1251 where we have used instead of 



10 Note that TJM11 allowed the ages to run over 15 Gyr in ex- 
trapolation, while here we impose 15 Gyr as a maximum age. It 
is also interesting to note that the SED-fit is a bit sensitive to 
the adopted IMF, and the Kroupa one gives slightly smaller stan- 
dard deviations, which makes sense. Standard deviations - after 
removal of the two outlying clusters NGC 5946 and 6544, are: age 
= 2.664; [Fe/H] = 0.226 for a Salpeter IMF; age = 2.51; [Fe/H] 
= 0.228 for a Kroupa IMF. 
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Best fit model: 7 Gyr, [Fe/H] = -1.3 
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Figure 22. GC SEDs from [Schiavon et al.l ||2005| . black), overplotted with their best-fitting, in terms of minimum 
X 2 , MILES-based template (red), for a Kroupa IMF, taking into account the HB morphology of each cluster. The 
age and metallicity parameters of the best fits are indicated and the residuals are given below each panel. 



the ages derived by De Angeli et al. (2005) , those by Marin- 
Franch et al. (2009) in which a correction for abundance ra- 
tio effects has been included. Though the basic conclusions 
remain unchanged, the CMD-ages shift towards older values. 
More quantitatively, the comparison with the CMD ages of 
De Angeli et al. (2005) is better for the SED-fit case, where 
the medians of the difference (CMD-age — integrated-age) 
are: -0.03± 1.8 Gyr for the SED-fit and -0.71 ± 2.83 Gyr 
for the TMJ. They become comparable, however, if we use 
the latest age determination by Marin- Franch et al. (2009), 
with corresponding values of 1.5 ± 2.2 Gyr 1.3 ± 2.76 Gyr, 
respectively. 

Far more striking is the comparison with the metallic- 
ity. Whereas in Figure [26] a remarkable agreement is evident 
between the Harris's (1996) metallicities and those derived 
with both the MILES-based SED-fitting and the Lick in- 
dex fitting, when we adopt the metallicity scale of Marin- 
Franch et al. (2009) (Figure US} the SED-fit metallicities are 
systematically lower. This may highlight that the SED-fit- 
derived metallicities mostly reflect [Fe/H] rather than the 



total metallicity [Z/H] that both Marin- Franch et al. (2009) 
and the TMJ models aim to trac e. This ambiguity in m etal- 
licity scales is amply discussed in lMaraston et all |2003l ) and 
iThomas et all (|201ll ). 

Finally, the values obtained with Lick indices directly 
calculated on the model SEDs without consideration of 
abundance ratio effects are off in any case (central panels 
in all figures). 

These figures highlight four important facts, namely: i) 
when using GCs as test for SSP models - the CMD-derived 
quantities need to be taken with care; ii) the analysis of ab- 
sorption line diagnostics should be carried out with models 
that account for element ratio effects otherwise the results 
are misleading; iii) selecting indices that - through models 
- are known not to be too sensitive to element ratio effects 
rectify the age and metallicity determination, which stresses 
again that the modelling of abundance ratio in the TMJ 
models is generally correct, and is required; iv) a good al- 
ternative appears to be the use of the full spectral fitting 
as the continuum shape helps getting ages and metallicity 
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Figure 24. Same as Fig[23] but for NGC 6544. Notice the large 
age-metallicity degeneracies, making this clusters one of the most 
poorly constrained. 
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Figure 26. As in Figure 1251 for the metallicity. On the y-scale 
we plot values from the Harris (1996 and updated versions) cat- 
alogue, which are held to reflect the [Fc/H] parameter. For the 
TMJ models we plot the total metallicity [Z/H]. 
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Figure 25. Comparison of ages obtained for the Schiavon (2005) 
GC sample as obtained via (panels from left to right): a) full 
SED fit on MILES-based SEDs (see Figure |22]l; b) Lick absorp- 
tion line indices calculated on the same MILES-based SEDs of 
case a), using: all available indices (filled symbols); a selection 
of indices the combination of which is little sensitive to element 
abundance ratios (TJM11) ; c) the TMJ Lick index models based 
on MILES with variable abundance ratios. On the y-scale we plot 
ages derived from CMD fit from De Angeli et al. (2005). 

correctly. Methods iii) and iv) cannot however release the 
abundance ratio parameter. 

We have shown in Section [575l that the empirically-based 
models are solar-scaled at high metallicity. In order to check 
the models in this regime, we fit the composite spectrum 
of the old open cluster M67 ijSchiavon et al.l[2004l ). which is 
believed to have a metallicity about or just slightly below 
solar, solar-scaled element abundance ratios, and a CMD fit- 
ting age of 3 to 5 Gyr depending on the adopted isochrones 
and fitting details (see Chaboyer et al. 1999 for a review and 
Schiavon et al. 2004 for full references). Figure [29] shows 
the best fit MILES-based and STELIB-based models for 
M67. The upper panels show the results for the minimum 
chi-square. STELIB-based models release a solution that is 



Figure 27. Same as Figure 1251 compared now with ages derived 
from CMD fit from Marin- Franch et al. (2009), where abundance 
ratio effects are explicitly considered. 

identical to the one from CMD-fitting. MILES-based models 
release a superb fit, which however gives the quite older age 
of 9 Gyr. The lower panels show the second ranked solution 
for STELIB-based and the 3rd ranked one for MILES-based 
models, which finally gets to the CMD solution. Note the 
risible difference in x 2 - 

Figure [3TJ1 shows the \ 2 distributions for age and metal- 
licity as obtained with STELIB-based and MILES-based 
libraries. Obviously, taking into account a 68% likelihood 
range instead of considering only the minimum x 2 would 
include the CMD solution for both type of models. We fur- 
ther verified that - using only solar metallicity MILES-based 
models - we obtain a best-fit age of 3 Gyr. Finally, we ver- 
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Figure 28. As in Figure 1271 for the metallicity published by 
Marin-Franch et al. (2009) which, after their correction for 
abundance-ratio effects is held to reflect the total metallicity 
[Z/H], 
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Figure 29. Model fits (red) to the composite spectrum of the solar metallicity old open cluster M67 (black, from Schiavon et al. 2004). 



ified that using the STELIB models tied onto the MILES 
temperature scale still releases the CMD solution, meaning 
that the solution is not driven by the near-IR slope. 

These results can be commented in several ways. One 
may drawn from here that the STELIB library better traces 
stars at solar metallicity, while MILES is a little too blue, 
which is why an older model is preferred. On the other hand, 
Vazdekis et al. (2010) claim to obtain the CMD solution as 
the best-fit. Here we note that the spectrum of M67 is not an 
observed integrated spectrum, but is a composite obtained 
by patching together individual stellar spectra via the pop- 
ulation synthesis code of Schiavon et al. (2004). Differences 
between their code and ours, both in terms of method and 
adopted isochrones, could couple just with the MILES li- 
brary in such a way as to explain why our first ranked solu- 
tion is not consistent with theirs. The Vazdekis et al. (2010) 
model is similar to the Schiavon's one (i.e. isochrone syn- 
thesis on Padova isochrones). This event, though unlikely, is 
not excluded. 

We refrain from drawing any conclusion on the libraries 



based on just one object, but we shall keep this fact in mind 
when fitting galaxies. 

What Figure [25] shows is that all fits are excellent and 
hardly distinguishable, which illustrates two aspects of spec- 
tral fitting for age/metallicity determination, namely that: i) 
the sole consideration of the minimum \ 2 is misleading (see 
Maraston et al. 2010 for a similar argument for the SED-fit 
of high-redshift galaxies) ; ii) the age-metallicity degeneracy 
is perfectly in place and the higher resolution of the models 
does not suffice to remove it. 



6.1.1 Effect of different data on parameter derivation 

The overlap b etwee n the iPuzia et al.l ([2002) and the 
ISchiavon et ail l|2005l) samples - consisting of eleven star 
clusters - allows a comparison between parameters derived 
using the same models but different observations. This is 
shown in Figures [31] and 1321 In most cases the outcome 
is the same, or at least a very similar, set of parameters. 
A few discrepancies can be seen, but considering that for 
these clusters there are noticeable differences in the shape 
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Figure 30. \ 2 distribution of age and metallicity obtained for 
M67 with MILES-based (lower panels) and STELIB-based (upper 
panels) models. 

of their spectra, and considering that even subtle differences 
can have an effect on the fitting, since age-metallicity degen- 
eracies are large and \ 2 distributions poorly constrained, the 
agreement is fairly good. The only exceptions are NGC 6218, 
whose SED shape in the P02 sample looks corrupted - us- 
ing the S05 spectrum we are able to recover parameters in 
good agreement with the literature values; the rebellious 
NGC 6553, for which the S05 version has been assigned 
a lower metallicity by 0.5 dex, while the P02 version has 
been given a higher metallicity by the same amount; and 
NGC 6626, where once again the metallicity derive d from 
the S 05 spectrum is in better agreement with the lHarrisI 
|l99rj ) value. 

6.2 Globular cluster B — V colours 

We compare GC colours with the models to check whether 
the new spectra help in removing existing discrepancies. No- 
tice that we limit ourselves to the B — V colour as the models 
show a clear trend, whereas for the other colours there is ei- 
ther no change, or not a clear one (as e.g. for the U — B) or 
existing data are not sufficient (as for the near-/i? colours). 
Figure [32] shows the B — V colours of MW GCs as a func- 
tion of metallicity, in a remake of the correspondent panel in 
Figure 21 of M05. Plotted are results for the MILES-based 
models (green) and for the M05 models (black) . Consistently 
with the discussion in Section 5.3, the B — V colours of the 
MILES-based models are blue r by 0.05 mag, a nd better fit 
the average locus of MW GCs. iMarastonl (|2005l ) anticipated 
this result indicating a similar offset to be applied to the 
M05 models and i t is neat to be a ble to find a reason for it. 

Furthermore, iPeacock et ail l|201ll) show that the q — 
r colours of M31 globular clusters are also much better 
matched by the new models because of the same effect, and 




A s e P02 ( G y r ) 

Figure 3 1. Comparison b etwe en ages derived for t he cross- 
sample of iPuzia et al. I ll2002l l and lSchiavon et al. I j2005l l globular 
clusters using SED-fitting with MILES-based models. 
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Figure 32. Same as Fig l31l but for metallicities. 



ISaglia et all (|2010h use the colour corrections and find a bet- 
ter match to the population of the old M31 Bulge. These 
results strengthen our conclusions regarding the mismatch 
between the Sloan colours of L RG galaxies and population 
models l|Maraston et al.ll2009bl . and Section 5.3). 
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Figure 33. B — V colours of M05 an d MILES-based models as a function of metallicity for 12 Gyr of age (and a Salpeter IMF) compared 
to Milky Way GCs data Jrlarrijl 996). This plot is identical to the correspondent panel of Figure 21 in M05. Note the improved agreement 
between models and GCs colours. 



7 SUMMARY AND CONCLUSIONS 

Stellar population models find ubiquitous applications in as- 
trophysics. The derivation of galaxy and star cluster prop- 
erties such as ages, metallicities, star formation histories, 
stellar masses via decoding of the integrated energy emis- 
sion is a crucial step in studies of galaxy formation and evo- 
lution and observational cosmology. Progress in technology 
now permits to obtain high spectral resolution for galaxies 
over a large range in masses and evolutionary status, i.e. 
age and level of on-going star formation, and for a relatively 
large range in redshift. Stellar population models have to be 
upgraded to meet the observational demand, and this was 
part of the aim for this work. 

The models have three main ingredients stemming from 
stellar evolution, namely the stellar energetics, the stellar 
atmospheric parameters and the individual stellar spectra, 
and also depend on the algorithm adopted in the code and 
on a variety of other assumptions (e.g. the IMF, the HB 
morphology, the mass- loss parameters, etc., see Maraston 
2005 for a wide discussion). The spectral resolution of the 
models is linked to the assumed individual stellar spectra. In 
this work we have varied the input stellar library - exploring 
a wide set of available empirical stellar libraries - keeping all 
other assumptions as in the M05 modelling. 

We replenish missing stellar types in the libraries with 
theor etical spectra from the MARCS |Gustafsson et all 
2008) theoretical library. We additionally calculate stel- 
lar population models fully based on the MARCS library 
which can be used to extend the wavelength baseline of the 
empirically-based models to longer wavelengths, all along 
to the AT-band and beyond. Similarly, we extend the mod- 
els into the UV by using our previously calculated high- 
resoluti on models based on the theoretical UVBLUE-Kurucz 
library (|Maraston et alj|2009bh . 

This is not the first work of this kind. Bruzual & Char- 
lot (2003) presented the first set of higher resolution mod- 
els exploiting the STELIB and the Pickles (1998) libraries, 



which have been the basis for many galaxy evolution studies. 
Vazdekis (199 9) calculated models using th e empir ical stel- 
lar library by I Jones! (119991 ). [Maraston et all (|2009al ) present 
very high-resolution SEDs up to ~ 4500 A using theoreti- 
cal Kurucz-type spectra. Vazdekis et al. (2010) provide their 
models based on the MILES library. 

The difference between this work and other similar in 
the literature is that - by exploring a wide range of empirical 
stellar libraries - we could also assess their effect on a stellar 
population model, which could be regarded as the second 
goal of this work. 

In general, the spectra of stars in common to the var- 
ious libraries present discrepancies in selected spectral re- 
gions, which are shown and discussed, most of which get 
smoothed out in the integrated stellar population models, 
which display - at the end - a fair level of homogeneity. 
However, some differences remain, which will be interesting 
to test with data of galaxies and star clusters. For exam- 
ple, MILES-based population models display lower flux at 
their near-IR edge - long-ward of ~ 6000A - with respect to 
population models based on both the empirical STELIB or 
Pickles libraries as well as to those based on the theoretical 
library MARCS. We could trace this effect as mainly due to 
the temperature scale adopted for RGB stars in MILES. 

As a relevant result, we confirm our earlier finding that 
the visual region of the spectrum of models based on em- 
pirical spectral libraries has lower flux (in a normalised 
sense) with res pect to the BaSeL-Ku rucz or the original 
Kurucz library (Ma raston et a l. 2009a). The only exception 
are the models based on the ELODIE library, that however 
have such a high resolution that flux calibration becomes 
more difficult. This effect helped to a better match of the 
g, r, i colours of Luminous Red Galaxies at redshift 0.4 and 
here we show that it also improves the comparison with the 
B - V vs [Z/H] relation for Milky Way GCs. This conclu- 
sion is in agreement with the finding of Peacock et al. (2011) 
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that the models of this paper better trace the g, r colours of 
globular clusters in M31. 

Most importantly, in this work we find that the the- 
oretical library MARCS agrees well with the empirical li- 
braries. This guarantees that the effect does not originate 
from the complicated procedure required to insert an em- 
pirical library in a population synthesis code. The adopted 
procedure matters, however. We have used the stellar pa- 
rameters provided by the libraries and scaled the empirical 
spectra using the theoretical flux in a well-defined region. 
Alternative approaches match the empirical spectra to the- 
oretical ones to insert the empirical libraries in the code (as 
done by Bruzual & Chariot 2003, for example). In this way, 
differences between empirical and theoretical spectra are di- 
luted. 

We also show that other recent models based on MILES 

- namely those by Vazdeki s et al. (2010) and by C onroy & 
Gunn (2010) - confirm the iMaraston et al.l (|2009eJ ) results, 
i.e. the impact of using empirical atmospheres on population 
models, even if details maybe different. Curiously, Conroy & 
Gunn (2010) came to the opposite conclusion. 

We also discuss extensively the Lick indices calculated 
directly on the integrated MILES-based SEDs and compare 
them with those calculated using the fitting function (FF) 
approach as well as with element ratio sensitive index models 
(Thomas, Maraston & Johansson 2011). In general, we find 
an excellent agreement between FF-based and SED-based 
index models. We also find a good agreement with the mod- 
els by TMJ, if the metallicity dependent chemical pattern of 
the Milky Way stars is properly taken into account in this 
comparison. SED-based models agree with solar-scaled TMJ 
models at high-metallicity and with [a/Fe]— enhanced mod- 
els at sub-solar metallicity, because stellar libraries carry 
into the models the spectral pattern imprinted of the ele- 
ment abundance ratios of the Milky Way field stars they are 
built upon. 

We tested the age and metallicity scales of the mod- 
els by fitting the observed SEDs of Milky Way GCs with 
ages and metallicities independently known from isochrone 
CMD fitting and resolved spectroscopy. We find overall good 
agreement when we perform full spectral fitting, suggest- 
ing that this an effective method to derive ages of stellar 
systems. This is because the spectral shape is age-sensitive 
even if details of spectral lines are not well-matched due 
to abundance ratio effects. Moreover, as we already noted 

- these models are - qualitatively - [a/Fe] -enhanced at sub- 
solar metallicities and solar-scaled at Zq and above, which 
also partly explains the good match to Milky Way GCs. For 
example, we could obtain excellent fits of the observed SED 
of the old open cluster M67 of the Milky Way disk, which 
has approximately a solar metallicity and 3 Gyr of age, and 
solar scaled abundance ratios. 

However, we find that - when we use - instead of the 
full spectral fitting - selected absorption lines without in- 
cluding a treatment of abundance-ratio effects - ages are 
underestimated and metallicities are overestimated. This is 
because absorption lines are very sensitive to element abun- 
dance ratio effects and peculiarities, which can be different 
between the stars composing the empirical libraries and ex- 
ternal stellar systems such as star clusters and galaxies. A 
better agreement is achieved when one uses a sub-sample 
of indices the combination of which is little sensitive to ele- 



ment abundance ratios (namely H« A , Mgb, Fe4383, Fe5335, 
Fe5406, see Thomas, Johansson & Maraston 2011 for de- 
tails). These facts should be considered when comparing 
these models to data. Using absorption line models including 
abundance ratio effects instead provide ages and metallici- 
ties in good agreement with the CMD-based ones (Thomas, 
Johansson, Maraston 2011) and the full spectral fitting per- 
formed with the SED models presented here. 

Finally, for M67 we noticed an interesting effect arising 
from the use of different libraries, namely that the mod- 
els based on STELIB releases as best-fit (minimum x 2 ) & 
solution which is identical to the one from CMD fitting, 
whereas the MILES-based models give an older age and a 
lower metallicity, even if the solution with parameters as 
from the CMD fitting has a \ 2 which is comparably good. 
This difference may originate from the MILES-based mod- 
els being slightly bluer than the STELIB-based ones, hence 
releasing older ages. Obviously, a result based on just one ob- 
ject is merely a speculation, but it will be interesting to check 
whether such an effect persists using galaxy data, which we 
shall pursue in a future paper. 

Models are available at www.maraston.eu/Mll. 
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APPENDIX A: REVISING THE NEAR-IR 
SLOPE OF MILES-BASED MODELS 

Longward 6300 A, the MILES-based models brandish lower 
flux levels with respect to the other models and to the 
MARCS-based models (cf. Figure [T5]). 

We traced this behaviour to have three origins: 

(i) The temperature calibration. For each iPicklesI (|l998l ) 
standard spectrum (we remind the reader that these are av- 
erages of several different stars of the same spectral type 
and luminosity class, and as such, should be less vulnerable 
to stochastic variations) we obtained the best-fitting MILES 
spectrum in terms of minimising the sum of squared residu- 
als. In general, it appears that the MILES temperatures are 
slightly underestimated with respect to the Pickles temper- 
atures. In particular, this holds true for Red Giant Branch 
bump and tip stars, whose impact on the integrated near- 
IR light of old stellar populations is significant. A spectral 
library with, on average, cooler temperatures, will have the 
effect that the integrated SED becomes bluer. This connects 
nicely to the discussion about the temperature scale in Sec- 
tions 13.51 and 3.6. Thus, by replacing the MILES tempera- 



tures with the Pickles ones for these stars, we are able to 
achieve higher flux levels red- ward of 6200 A. 

(ii) Gaps in stellar parameter space. By shifting the 
coolest giants to slightly higher temperatures, there are no 
longer any cool enough stars to cover the parameter space 
required by the adopted stellar evolution prescription for 
the tip of the RGB and E-AGB. By adding sufficiently cool 
stars, for example from the Pickles library, part of the flux 
drop can be cured. 

(iii) Short wavelength range systematic offset. Besides the 
general lower flux long-ward 6300 A we spotted another 
near-IR feature on MILES spectra, namely a small dip in the 
flux levels between 6500 and 7000 A. This feature emerged 
when we compared MILES spectra - smoothed to the ap- 
propriate resolution - with Pickles spectra. We found them 
to agree perfectly over the entire wavelength range, except 
in that short interval. This discrepancy can certainly not be 
ascribed to any difference in temperature, gravity, etc., but 
must be an intrinsic difference between Pickles and MILES 
spectra. Since the Pickles spectra are averages of several ob- 
servations from numerous different sources, while the MILES 
spectra have all been obtained using the same instrumental 
setup, the most likely explanation is that this is an artefact 
originating from the MILES instrumental setup. Alterna- 
tively, this could be due to the adopted procedure for cor- 
recting against telluric absorption. Note that we spotted this 
mini-dip only on some MILES spectra as - in order to notice 
it - one needs to have an otherwise identical spectrum at 
moderate resolution, such as a Pickles spectrum. The dip is 
probably present in all MILES spectra - we argue it is a sys- 
tematic - but we could not check all spectra in this respect 
for the above named reasons. 

Fig. IA 1 | shows such best-fitting MILES spectra to four 
standard IPicklesI { 1998! ) spectra, accompanied by their re- 
spective flux ratios as a function of wavelength. MILES spec- 
tra have been smoothed and rebinned to the Pickles approx- 
imate resolution of R«500, and 5A sampling. It is obvious 
that, even when the two spectral energy distributions agree 
on general very well, there is a systematic deviation rewards 
~ 6500 A. To counterbalance this, we have chosen 12 of the 
best-fitting MILES spectra (6 dwarfs and 6 giants), fitted 
second-order polynomials to the flux ratios in the aforemen- 
tioned spectral range for each of the 12 components, and ap- 
plied the straight average of these fitting curves (see Fig |A"2"1l 
to all the MILES stars. 



Based on the above corrections, for the solar metallicity 
case we provide two versions of our MILES-based models, 
original and revised. Figure A3 shows how the correspondent 
spectra compare, for the illustrative age of 12 Gyr (all other 
ages behave the same way). One notice from the residual plot 
that - since the most affected stars are cool giants - the near- 
IR portion of the spectrum is the only one that is affected 
appreciably by our revision. We recall that the same exercise 
at different metallicities cannot be made as it is based on 
Pickles spectra and they exist mostly at solar metallicity. 
It will be interesting to understand the physical meaning 
of such revision - if any -, which can be accomplished by 
comparing these models with stellar, star cluster and galaxy 
data. 



High-resolution stellar population models 
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Figure Al. Best-fitting MILES spectra to four standard spectra from Pickles (1998), and their flux ratios as a 
function of wavelength, defined as AF\ = F\ pickles/ F\ MILES (smaller panels). Flux ratios are scattered around 1, 
except redwards of ~6500 A, where a systematic flux deficit in the MILES spectra can be observed. For this exercise 
the MILES spectra have been smoothed and rebinned to the Pickles resolution and sampling. 
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Figure A2. Second-order polynomial fits to the flux ratio plots 
of best-fit ting MILES sp ectra to 6 dwarf and 6 giant standard 
spectra of IPickled 1)19981 1. in the wavelength range 6450-7050 A 
{dashed black lines). The straight average applied to all MILES 
stars is also shown {red solid line). 




Figure A3. Comparison between MILES-based and MILES- 
revised based solar metallicity SSPs for the illustrative age of 12 
Gyr The difference in flux (shown in the lower plot as F rev - ma ^ — F) 
induced by our revision is very mild all over the spectrum and 
mostly affect the near-IR part. 



